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A deep Chandra observation of the poor cluster AWM 4-1. 
Properties of the central radio galaxy and its effects on the 
intracluster medium 
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ABSTRACT 

Using observations from the Chandra X-ray Observatory and Giant Metrewave Radio Tele- 
scope, we examine the interaction between the intracluster medium and central radio source in 
the poor cluster AWM 4. In the Chandra observation a small cool core or galactic corona is re- 
solved coincident with the radio core. This corona is capable of fuelling the active nucleus, but 
must be inefficiently heated by jet interactions or conduction, possibly precluding a feedback 
relationship between the radio source and cluster. A lack of clearly detected X-ray cavities 
suggests that the radio lobes are only partially filled by relativistic plasma. We estimate a 
filling factor of 0=0.21 (3(t upper limit < 0.42) for the better constrained east lobe. We 
consider the particle population in the jets and lobes, and find that the standard equipartition 
assumptions predict pressures and ages which agree poorly with X-ray estimates. Including 
an electron population extending to low Lorentz factors either reduces (7„im = 100) or re- 
moves (7mm = 10) the pressure imbalance between the lobes and their environment. Pressure 
balance can also be achieved by entrainment of thermal gas, probably in the first few kilopar- 
secs of the radio jets. We estimate the mechanical power output of the radio galaxy, and find 
it to be marginally capable of balancing radiative cooling. 

Key words: galaxies: clusters: general — galaxies: clusters: individual (AWM 4) — inter- 
galactic medium — galaxies: active — cooling flows — X-rays: galaxies: clusters 



1 INTRODUCTION 

X-ray observations of clusters and groups of galaxies over the 
last decade have led to a significant revision of our models of the 
intergalactic medium in these systems. The Chandra and XMM- 
Newton observatories have provided strong evidence that despite 
cooling times being sig nificantly shorter than the Hubble time (e.g., 
ISanderson et alj 20061), relat i vely little gas act ually cools below 



4 keV l lPeterson et al.l2003l : lKaastra et alj2004l) . It is now widely 
accepted that excessive cooling is in many systems prevented by 
a feedback mechanism in which the AGN of the central domi- 
nant galaxy, fuelled by cooling intra-cl uster gas, can reheat the 
gas through a variety of me chanisms (e.g.. |Peterson & Fabiai]|20()q : 
iMcNamara & N ulsen' 2007. and references therein). 

X-ray and radio images provide numerous examples of inter- 
actions between AGN and the surrounding intra-cluster medium 
(ICM). Deep multi wavelength observations of the brightest nearby 
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clusters have revealed complex structures associated with the ra- 
dio jets and lobes, including shocks, sound waves, individual and 
linke d chains of ca vities, uplifte d material and cooli n g filaments 
(e.g. [Fabian et al.l Eo OS. 2006; iForman et all |2007| ; IWise et al.l 
l2007l : iBlanton et aljl2009.) . Much of this work has concentrated 
on the cavities in the ICM which radio lobes produce as they in- 
flate. The enthalpy of the cavities can be used as a measure of the 
mechanical power output of the radio jets, and has been shown to 
be sufficient to pre vent or greatly reduce cooling in many systems 
teirzan et alj|2004 l). provided the energy can be transfered into the 
intracluster medium and distributed relatively isotropically. 

The disturbed structures produced by AGN jet/ICM interac- 
tions are relatively short-lived, and increasingly difficult to detect 
as they age. The radio lobes which inflate cavities fade rapidly once 
the AGN outburst ceases, and the X-ray cavities, which are de- 
tected by contrast with their surroundings, become less visible once 
they move beyond the dense group or cluster core. It is therefore 
considerably more difficult to study older AGN outbursts. How- 
ever, since the mechanism by which cavities heat their surround- 
ings is still a matter of debate, it is desirable to observe systems 
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with as wide a range of ages as possible, so as to understand clearly 
the interaction between the radio lobes and their environment. 

Most observations of radio galaxies to date have been based on 
observations at frequencies > 1 GHz. Radio lobes may be studied 
over a wider range of timescales by observations at lower radio 
frequencies, which probe lower energy electrons less affected by 
spectral aging. Deep X-ray imaging is needed to complement such 
observations, and in this paper we discuss one example where this 
combination is available, the poor cluster AWM 4. 

A previous XMM-Newton observation found AWM 4 to be ap- 
proximately isotherma l to a radius of ^^150 k pc, with no evidence 
of a central cool core jO'SullivanetalfcOOl hereafter referred to 
as OS05). Comparison with MKW 4, a cluster of similar tempera- 
ture and galaxy population, but which hosts a large cool core and 
lacks a central radio source (O'Sullivan et al. 2003), leads to the 
suggestion that AWM 4 had been strongly heated by its central ra- 
dio galaxy, 4C-I-24.36. However, examination of the XMM-Newton 
data showed no evidence of cavities or shocks, and no high reso- 
lution images showing the lobes were available in the literature or 
radio archives. The existing VLA 1.4 GHz data were interpreted 
as evidence against the presence of lobes of su fficient volume to 
be re sponsible for reheating a large cool core jGastaldello et al.l 
l2008l) . The ICM and galaxy distribution both appear relaxed with 
no significant substructure, wi th a strong concentratio n of early- 
type galaxies toward the core dKoranvi & Geiie3l2002h . A cluster 
merger, which could also have heated the ICM, therefore appears 
unlikely. The central elliptical, NGC 6051, shows no signs of re- 
cent interactions (Schombert 1987) and is considerably more lumi- 
nous than its neighbours, with a difference in magnitude above the 
second-ranked galaxy of i\fi2=1.6 (SDSS g-band). 

An analysis of deep GMRT radio observations at 235, 327 and 
610 MHz was presented in (Giacintucci et al. 2008, hereafter re- 
ferred to as GVM08). These data provided much new information 
about 4C-I-24.36, revealing radio emission from the jets and lobes 
extending ^^^75 kpc from the AGN. The source was shown to be a 
wide-angle-tail radio galaxy with inner jets oriented close to the 
plane of the sky, probably moving southward with a velocity of 
<120 km s~^. From modelling of the progressive steepening of the 
spectral index a (defined as 5 oc where 5* is flux and v fre- 
quency) along the jets the radiative age of the electron population 
was estimated as 160-170 Myr. 

In this paper, we use a new Chandra ACIS-S observation of 
AWM 4, in combination with the GMRT and archival VLA obser- 
vations, to study the structure of the ICM and the interaction of the 
AGN, radio jets and lobes with the surrounding hot gas. The gen- 
eral properties of the cluster are summarised in Table[T] along with 
the position, distances and angular scale of the system. Throughout 
the paper we assume Ho=70, Q.m = 0.3, and 57a = 0.7. Un- 
certainties are generally quoted at the la level, except in the case 
of X-ray spectral fitting, where 90 percent uncertainties were esti- 
mated. Section|2]describes the observation and data reduction, and 
Sections |3] and |4] our imaging and spectral analysis. In Section |5] 
we examine the properties of the gas in the core of NGC 6051, 
immediately surrounding the AGN, and in Section |6T| we discuss 
the interaction between the radio jets and ICM, and place limits on 
the timescale of the outburst and the particle content of the radio 
lobes. We discuss our results in Section|7]and list our conclusions 
in Section[8] 



Table 1. General properties of the AWM 4 system 
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S235MHZ (mJy) 


2750 




logPiAGHz (W Hz-1) 


24.14 



2 OBSERVATIONS AND DATA REDUCTION 

AWM 4 was observed by the Chandra ACIS instrument during 
Cycle 9 on 2008 May 18-19 (ObsID 9423) for ~80 ks. A s\im- 
mary of the Chandra m ission and instrumentation can be found in 
Weiss kopf etakl ( l2002h . The S3 CCD was placed at the focus of the 
telescope and the instrument operated in very faint mode, to take 
advantage of the superior cosmic ray rejection. We have reduced 
the data from the pointing using CIAO 4.1.2 a nd CALDB 4.1.2 
following techniques similar to those described in IO' Sullivan et all 
12007) and the Chandra analysis thread|3- The level 1 events files 
were reprocessed using the very faint mode filtering, bad pixels 
and events with ASCA grades 1, 5 and 7 were removed, and the 
cosmic ray afterglow correction was applied. The data were cor- 
rected to the appropriate gain map, the standard time-dependent 
gain and charge-transfer inefficiency (CTI) corrections were made, 
and a background light curve was produced. The observation did 
not suffer from significant background flaring, and the final cleaned 
exposure time was 74.5 ks. While data from the entire detector were 
examined, for the purposes of this study we generally only make 
use of the S3 CCD, as the radio source and the cluster core fall on 
that chip. 

Identification of point sources on S3 was performed using the 
CIAO task WAVDETECT, with a detection threshold of 10"", chosen 
to ensure that the task detects ^ 1 false source in the field, work- 
ing from a 0.3-7.0 keV image and exposure map. Source ellipses 
were generated with axes of length 4 times the standard deviation of 
each source distribution. These were then used to exclude sources 
from most spectral fits. A source was detected coincident with the 
peak of the diffuse X-ray emission; this was considered a poten- 
tially false detection and ignored, though we did later test for the 
presence of a central X-ray point source. 

Spectra were extracted using the SPECEXTRACT task. Spectral 
fitting was performed in XSPEC 11.3.2ag. Abundances wer e mea - 
sured relative to the abundance ratios of lGrevesse & Sauva]| ( ll998h . 
A galactic hydrogen column of 0.05 x lO^'^ cm~^ and a redshift of 
0.0318 were assumed in all fits (except those whose purpose was 
to test the effect of varying hydrogen column), and 90% errors are 
reported for all fitted values. Spectra were grouped to 20 counts per 
bin, and counts at energies above 7 keV were ignored during fitting. 

Background spectra were drawn from the standard set of CTI- 
corrected ACIS blank sky background events files in the Chandra 
CALDB. The exposure time of each background events file was 
altered to produce the same 9.5-12.0 keV count rate as that in the 
target observation. The same very faint mode background screen- 
ing was applied to the background data sets. Comparison of source 
and background spectra suggested a slight excess of soft emission 

1 http://asc.harvard.edu/ciao/threads/index.html 
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in the background datasets, mainly below 0.5 keV. This is not unex- 
pected, as the soft X-ray background arises largely from hot gas in 
the galaxy, and from coronal emission associated with solar wind 
interactions, an d thus is bo th spatially and temp orally variable (e.g., 
iKuntz & Snowd en 2000; Snowden et al .120041) . There are also indi- 
cations that the spectral shape of the background has changed since 
the creation of the blank-sky background files (c.f. the ACIS back- 
ground cookbool|3), which could contribute to the disagreement at 
low energies. Experimenting with fitting spectra in different energy 
bands, we found that ignoring energies below 0.7 keV produces re- 
sults which are consistent with those derived from XMM-Newton in 
OS05. If the absorbing hydrogen column is allowed to vary in these 
fits it still tends to produce best fit values in excess of the galactic 
value, but typically consistent with it at the 90% uncertainty level, 
and with the hydrogen column derived from the XMM-Newton fits. 
Temperature, abundance and normalisation of the models are not 
significantly affected if the hydrogen column is allowed to vary, 
and our results are therefore independent of its value. 



3 IMAGING ANALYSIS 

An initial examination of exposure-corrected images of AWM 4 
agrees with the general conclusions of the XMM-Newton anal- 
ysis of OS05. The cluster appears to be relaxed, with a fairly 
smooth elliptical surface brightness distribution roughly centred on 
NGC 605 1 . Figure[T]shows a soft-band image of the cluster with the 
D25 ellipse of NGC 6051 and GMRT 610 MHz contours overlaid. 
There are no surface brightness discontinuities or "fronts" which 
would indicate the presence of large-scale gas motions or shocks 
in the cluster halo. However, a number of features which were not 
observed in the XMM images are visible in the Chandra data. The 
most obvious of them is a bright central peak, closely aligned with 
the central cD and the radio core detected at 4.9 GHz. This is only 
a few pixels in diameter, and so could not have been resolved by 
XMM. We will discuss this in more detail in Section [3T| 

Weak features are visible along the jets, particularly in the east 
lobe and the knots ~20 kpc from the radio core. To examine these, 
we created unsharp masked images smoothed at a range of scales. 
Figure |2] shows one example. The strongest feature in the field is a 
surface brightness deficit in the east radio lobe. The deficit is not 
well-correlated with the radio structure, instead appearing as two 
small "holes" in the X-ray emission toward the end of the lobe. No 
surface brightness structures are seen in the western lobe, though 
there is some indication of a broad opening or bay in the surface 
brightness on the west side of the cluster, into which the jet flows. 

Using the CIAO SHERPA application we performed a two- 
dimensional surface brightness fit to the S3 0.7-3 keV image and 
subtracted the model image to search for residual features. Test- 
ing showed that PSF convolution did not significantly affect the 
fit, and we therefore used an unconvolved model. The background 
contribution was estimated from the scaled blank-sky dataset, and 
we corrected for vignetting and other effects using a monoener- 
getic 1.05 keV exposure map, the energy chosen to match the mean 
photon energy of the data. A reasonably accurate fit was produced 
with two /3-models, an extended elliptical component describing 
the large-scale emission and a compact component with negligible 
ellipticity describing the core surface brightness peak. The best fit 
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Figure 1. 0.3-2 keV exposure-corrected Chandra image of AWM 4, 
smoothed with a Gaussian of width 5 pixels. The dashed blue line marks 
the R-hrnid D25 ellipse for the cD galaxy, NGC 605 1 . Contours (coloured 
black or white for clarity) represent the GMRT 610 MHz radio intensity, 
taken from GVM08. and are spaced by factors of 2, starting at 0.2 mJy b^ ^ . 
The restoring beam was 5" X 4" . 



Table 2. Parameters of 2-d surface brightness model used to create residual 
map 



Component 


Parameter 


Value (1(T error) 


1 
2 


rcore ( ) 

/3 

T^core ( ) 


14+0-" 
(,400+0.018 

39-90 




/3 


0.403±0.002 




ellip. 


0.23I±0.007 




p.a. (°) 


80.2±0.8 



Position angle (p.a.) represents the angle of the major axis north of due 
west. 

parameters are given in Table[2]and the residual image is shown in 
Figure |2] 

Unfortunately for our purposes, the strongest features in the 
resulting residual map are associated with a general inequality in 
surface brightness between the north and south quadrants of the 
cluster. This is visible in the residual image as a strong deficit in sur- 
face brightness to the south of the radio jet. AWM 4 is imperfectly 
described by a simple elliptical model, owing to the presence of 
additional surface brightness to the north of the core. This suggests 
that NGC 6051 is located somewhat to the south of the centroid of 
the large-scale X-ray halo. GVM08 determined that NGC 605 1 is 
probably moving south, and this may explain the offset between the 
galaxy and ICM centroid. There is also an apparent arc of brighter 
emission to the north of the core, whose origin is unclear. 

The residual image does reveal some indications of structures 
associated with the jets and lobes, most notably X-ray bright re- 
gions along the boundaries of the jets and possibly the inner part of 
the eastern lobe. The broad opening or "bay" in the surface bright- 
ness on the western side of the cluster also produces some features, 
though the opening angle is considerably wider than that of the jet 
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Figure 2. Upper panel: Chandra 1-3 keV unsharp mask image, based on 
images bimied by factor 2 and smoothed with 3 and 30 pixel Gaussians. 
Lower panel: 0.7-3 keV residual image, after subtraction of best-fitting sur- 
face brightness model, smoothed with a 7 pixel Gaussian. 610 MHz radio 
contours are overlaid on both images (see Fig.[T]for details). Regions dis- 
cussed in the text are marked; ellipses and circles mark regions used for 
significance calculations, while dashed lines indicate possible structures in 
the X-ray emission. 



and lobe. All of these features are relatively weak and identified in 
part because of apparent correlations with the jets. 

To test the statistical significance of the features, we compared 
the structures with regions at similar radii, using the surface bright- 
ness models to define approximate isophotal ellipses. An elliptical 
annulus of width 40" covers the surface brightness deficits in the 
eastern lobe and the western bay, running through the base of the 
western lobe and the centre of the eastern lobe. We measured the 
1-3 keV exposure corrected surface brightness (with point sources 
removed) in an azimuthal profile around the ellipse, using sections 
with angular width 9°or 18°. Figure[3]shows the results. 

The two strongest negative features in the azimuthal profile 
occur within the boundaries of the radio lobes, and correlate with 
the features seen in the images. The deficit in the east lobe is 
stronger but smaller (corresponding to the "holes") while the west- 
em deficit is broader but less deep, and offset to the north side of 
the lobe. 

We estimate the significance of the deficits in the lobes by cal- 
culating the mean number of counts in the regions of the azimuthal 
profiles. Excluding the regions within the lobes (as marked on Fig- 
ure O the mean of the 9° regions (or 18°), exposure corrected and 
normalised to a fixed area of 4600 (2300) pixels, is 640.9±25.3 
(320.5±17.9) counts. The lowest value in the east lobe is 547.1 
(257.7) counts, giving a 3.7(7 (3.5(7) significance. This suggests that 
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Figure 3. 1-3 keV Azimuthal surface brightness profile around AWM 4 at 
the radius of the eastern lobe. Black points show 18° sectors, grey points 9° 
sectors. The solid line marks the mean surface brightness excluding the 3 
lowest black (6 lowest grey) points, dashed lines show the Icr uncertainties 
on the mean using each set. Dotted lines mark the angle within which the 
radio lobes are detected at 235 MHz. 



the "holes" in the east lobe represent a significant surface bright- 
ness depression, probably indicating the presence of a cavity in the 
ICM. 

The lowest points in the western lobe are 2.9(7 (2.3(7) signif- 
icant, and given the poorer alignment with the radio lobe it seems 
possible that they represent a more general deficit on this side of the 
cluster, rather than a cavity associated with the western lobe. This 
is supported by a simple analysis of the X-ray surface brightness 
in an ellipse chosen to match the lobe as a whole, compared with 
its surroundings. This shows a deficit only at Icr significance. Sim- 
ilar analysis for the east and west jet knots yields low significances. 
The highest, for the west knot, is 2.5a. Deeper observations would 
be needed to determine the nature of these apparent structures. 

The lack of clear, significant structures in the cluster is in it- 
self interesting. If the radio lobes have excavated cavities in the 
ICM, we would expect a surface brightness deficit on their line of 
sight. We can estimate the volume of the lobes from their projected 
size. We assume them to be oblate ellipsoids with major/minor axes 
21 X 14.5 kpc (West) and 22.3 x 12.8 kpc (East) at radii of 68.8 and 
64.7 kpc respectively. From the surface brightness models we ex- 
pect to see 1 120 counts (1-3 keV) along the line of sight of the west 
lobe and 1270 from the east lobe if there are no cavities. Based on 
the density profile of the ICM (see Section|4ll, we estimate that we 
would expect deficits of 184 counts and 25 1 counts for the west and 
east lobes respectively if the lobes are empty of ICM plasma. We 
would therefore expect to detect cavities at reasonable significance 
(4-5(7). The fact that we find a significant decrement only in a small 
region of the east lobe raises questions about the lobe contents. 



© 2010 RAS, MNRAS 000,[TlfT8l 



A deep Chandra observation ofAWM 4 5 




Figure 4. Unbinned Chandra images of the core of AWM 4 in four energy 
bands, with VLA 4.9 GHz contours overlaid. The energy bands are (up- 
per left) 0.3-1 keV, (upper right) 1-3 keV, (lower left) 3-5 keV and (lower 
right) 5-7 keV. The contour levels are spaced by factors of two starting at 
0.15 mjy b-i and are taken from an image with half-power beam width 
(HPBW) 0.7"x0.4". The images are ~12 kpc or ~19" wide. 
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3.1 The cluster core and AGN 

The observed central surface brightness peak could arise from sev- 
eral causes, including X-ray emission from the AGN, the presence 
of a small-scale cooling region or galactic corona, and emission as- 
sociated with the stellar population of NGC 6051. Surface bright- 
ness fitting to the 0.7-3 keV image confirms that the emission is 
extended and is poorly described by a point source model, indi- 
cating that an AGN cannot be the only source present. Figure |4] 
shows images in four energy bands, 0.3-1 keV, 1-3 keV, 3-5 keV 
and 5-7 keV. The central source is clearly extended up to 3 keV, but 
examination of the 3-5 keV image shows four counts detected at 
the position of the radio core, and no counts in the 5-7 keV band. 
Contributions from thermal plasma emission should be negligible 
in such a small region in the harder bands, so the 3-5 keV counts 
may arise from either AGN or X-ray binary emission. 

Comparison with the 4.9 GHz contours (derived from VLA 
archival data, project AK0360), indicates that the extent of the sur- 
face brightness peak is similar to the distance at which the jets are 
first detected, particularly for the western jet. The lack of detectable 
jets inside this region could indicate that they have switched off. 
We consider this unlikely, since the spectral index in the innermost 
region where the jets are detected is a ~ 0.5, consistent with a 
currently active radio galaxy rather than a dying source (GVM08). 
We instead assume that the jets are not detected closer to the core 
because they are highly collimated and therefore not resolved. This 
suggests that the jets lose collimation and 'flare' to a greater width 
close to the edge of the X-ray surface brightness feature. 

We can estimate the emission likely from unresolved point 
sources in NGC 605 1 based on known X-ray /optical scaling rela- 
tions. Unfortunately, the scale of the surface brightness peak is only 
a few arcseconds, and no suitable high-resolution optical data are 
available. We therefore extract a ATs-band surface brightness pro- 
file from the 2-Micron All-Sky Survey (2MASS, ISkrutskie etal] 
[2OO6) image of NGC 6051. This image has pixels 1" in size, but 
is derived from a combination of 2"-pixel images obtained by the 
survey telescope. The central /^s-band surface brightness profile is 
effectively smoothed by this process and may be broader and less 
peaked than the true profile. 



Figure 5. 0.5-2 keV surface brightness profile of the central 20" of AWM 4 
(errorbars indicate Icr uncertainties), compared to estimates of the contri- 
bution to surface brightness from point sources in NGC 605 1 and the intra- 
cluster medium (ICM). Surface brightness profiles for low mass X-ray bi- 
naries and for stellar sources (cataclysmic variables and coronally active 
binaries) are estimated by scal ing the 2MASS A"s-ban d sur face brightness 
profil e using the relations of iKim & Fabbiandl2004^ and jSazonov et alj 
I2OO6I) . The ICM contribution is determined from a /3-model fit to the emis- 
sion between 30 and 200". 1 ct uncertainties on the sum of these profiles are 
estimated from the uncertainties on the scaling relations and the normalisa- 
tion of the surface brightness fit to the ICM. 

We estimate the X-ray flux from low mass X-r ay binaries 
(LMXBs) using the relation of lKim & Fabbiand ( |2004|) . The small 
number of 3-5 keV counts observed at the position of the radio 
core are consistent with the expected LMXB emission in this area, 
within the (large) uncertainties. We place a 90% upper limit on the 
AGN luminosity of Lagn = 1.1 x 10*" erg s"^(0.7-7 keV) by as- 
suming a power law spectrum with r=1.7, and subtracting our best 
estimate of the flux from low-mass X-ray binaries. 

To calculate the expected point source contribution in the area 
of the extended soft-band core, we must also include emission from 
stellar sources such as cataclysmic variables and coronally active 
binaries. This contribu tion can be pre dicted using estimates from 
the Milky Way bulge Isazonover^ 2006). We work in the 0.5- 
2 keV band for which this relation is defined, scaling the LMXB 
contribution to match. The stellar sources are assumed to have a 
0.5 keV, solar abundance thermal plasma spectrum. We estimate 
the contribution from the large-scale gaseous halo of AWM 4 by 
fitting a /?-model to the 0.5-2 keV surface brightness profile be- 
tween 30" and 200". This produces an almost flat contribution to 
surface brightness in the core. Figure [5] shows these predicted sur- 
face brightness profiles compared to the measured surface bright- 
ness. 

At radii greater than 6", the sum of the predicted profiles is a 
reasonable match to the measured surface brightness profile. Inside 
this radius the agreement is less satisfactory; the measured profile 
is brighter than expected by a factor ~3 at the centre and appears to 
be more sharply peaked. The overestimation of flux at 2-6", while 
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statistically acceptable, is probably a result of the poorer resolution 
of the ifs-band profile. It seems unlikely that the underestimation 
of the peak flux could arise from this source, and we therefore con- 
clude that an additional extended spectrally soft surface brightness 
component is present. 



4 SPECTRAL ANALYSIS 

To establish the general radial dependence of gas properties, we 
extracted spectra from the S3 CCD in circular annular regions cho- 
sen to produce spectra with a signal-to-noise ratio greater than 100 
in a 0.7-7.0 keV band. A central bin of radius 6 pixels (~3") was 
added to allow us to search for differences in the central surface 
brightness peak. The annuli were centred on the peak of the X-ray 
emission, which corresponds to the optical centre of NGC 605 1 and 
the position of the cor e of 4C+24.36. The spectra were fitted using 
the WABS absor ption i Morrison & McCammonlll983l) and APEC 
s Is 



Radius (kpc) 

10 



plasma models dSmith etal.ll200lh . and since gas properties are 
expected to vary with true three-dimensional radius, we used the 
XSPEC PROJCT model to deproject the data. Beyond ~200" the 
annuli no longer fall entirely on the S3 chip. The lost area was cor- 
rected for in the PROJCT model, using the angular coverage param- 
eters for each annulus. The outermost annulus in each fit contains 
less than the required number of counts and is not deprojected; for 
this reason values in the outermost bin should only be considered 
indicative. The radio lobes only extend to ~145", so we do not 
expect edge effects to influence our results. 

A second set of annular spectra were also extracted, with the 
radii of the bins in the range 3-136" chosen to be closely com- 
parable with regions used in examining the physical properties of 
the radio jets and lobes. These spectra are treated identically to the 
S/N=100 profile, and produce very similar results. Comparison of 
the two profiles shows that outlying data points, such as the temper- 
ature in bin 5 of the S/N=100 profile, and abundance in bin 2 of the 
alternate profile, are not seen replicated when the radii of the spec- 
tral regions are changed. This suggests that they do not represent 
the true temperature or abundance at those radii; the underlying 
gas properties are likely to vary smoothly with radius. 

Figure |6] shows the deprojected temperature and abundance 
profiles for this deprojection. Comparison with the XMM tempera- 
ture and density profiles shows close agreement, confirming the ap- 
proximate isothermality and smooth decline of density with radius 
over the field of view. The exception is the central bin, which has 
a temperature of ~ 1.2 keV compared to the ~2.6 keV surrounding 
emission and a mean density approximately an order of magnitude 
greater than the ICM. This suggests that the central surface bright- 
ness peak corresponds to a small-scale cool core or galactic corona. 

An apparent difference in abundance profiles between Chan- 
dra and XMM arises fro m the different abundance tables used in 
the anal ysis. We use thelGreyesse & Sauvall ( [1993) tables, while 
the older 'Anders & Grevessd 1 19891) tables were used for the XMM 
analysis (OS05). Refitting our models using the older table pro- 
duces results in good agreement with the XMM profile. There is 
some indication of a central abundance plateau in AWM 4, extend- 
ing out to ^^75" ('^SO kpc). Outside this radius the abundance de- 
clines. 

The central bin has a very low abundance, possibly owing to 
our use of a single temperature plasma model for a spectrum which 
must have multiple emission components. We therefore repeat our 
deprojection analysis, adding two more components to the spectral 
model of the central bin, a 0.5 keV solar abundance APEC model 




1 10 100 

Radius (arcsec) 

Figure 6. Deprojected temperature, abundance and density profiles for 
AWM 4. Black diamonds indicate fits to Chandra spectra from annuli with 
signal-to-noise ratio >1()0, crosses to a modified set of annuli selected for 
comparison with the radio data. Diamond points ai'e omitted from the den- 
sity plot for clarity, but agree with the crosses within uncertainties. Grey 
rectangles show the XMM deprojected temperature and abundance profiles 
and the soUd line shows the density profile derived from the XMM data 
(with 1(T uncertainties indicated by the grey shaded region), all from OS05. 



representing stellar sources and a r=1.7 power law for the LMXBs. 
Normalisations were determined from the estimated surface bright- 
ness profiles described in Section lSTI These additional components 
affect the parameters of the original plasma model only slightly, 
and we find a temperature of fc_gr=1.00^Q 'j7 keV and abundance 
of O.IOIq '^Q Zq. This is still very low, and it seems likely we are 
still underestimating the abundance, possibly owing to the presence 
of a range of temperatures within the central bin (The Fe-bias ef- 
fect, .Buote & Fabian 1998.) . 



5 GAS PROPERTIES OF THE CORE 

The mean density in the core, assuming that it fills a sphere of ra- 
dius 1875 pc (-^3"), equ ivalent to our spe ctral extraction region, 
is ne=0.10lto 035 cm~^. ISunetalJ ( |2007|) examined a sample of 
small cool cores contained within galaxies, which they argued were 
formed largely from material produced by stellar mass loss, dis- 
tinct from the surrounding hotter gas associated with the larger 
group or cluster. They suggest that these galactic coronae should 
be in pressure equilibrium with the surrounding ICM. Estimat- 
ing the required pressure, we find that a density of 0.291 cm^^ is 
needed. If the background subtracted flux in the 3" aperture arises 
only from a corona with this density and the measured tempera- 
ture, it must occupy a smaller volume, a sphere of radius 1320 pc 
(^^2"). This is in reasonable agreement with our surface brightness 
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modelling. The 0.5-2 keV and 1.4 GHz luminosities, Lx, 0.5-2 = 
1.76 X 10''° erg s"^ and Lr,i.4 ghz = 1.31 x 10^* W Hz-\ are 
in the respective ranges found for the population of coronae studied 
by (Sun 2009). 

The isobaric cooling time of the gas in the cool core can be 
approximated as 

bknTneVfie 



tcool 



(1) 



where V is the volume of the gas, Lx,boi is the bolometric 
X-ray luminosity, fcsT is deprojected temperature, the depro- 
jected electron number density, and and /^e are the mean molec- 
ular weight (0.593) and mean mass per electron (1.167) of the gas. 
We find a cooling time of 302j;gg" Myr, as compared to tt,ooi=l-5- 
3 Gyr for the ~2.6 keV gas immediately outside the corona. 

Gas cooling from the corona may be the main source of mate- 
rial fuelling the active nucleus. We therefore consider the physical 
parameters of this cooling to determine whether the corona is a re- 
alistic source of fuel over the lifespan of the outburst. The mass of 
gas cooling out of the X-ray phase within the corona can be esti- 
mated based on the assumption of continuous isobaric cooling as 

2^,mpLx,boi 



Mcool 



5kBT 



(2) 



where nip is the proton mass. We find a deposition rate of 
Mi.ooi=0.067 M0yr~ . This is a sufficient rate to support the ob- 
served AGN activity; taking our largest estimate of the mechanical 
energy output of the AGN, 1.63±0.02 x 10^^ erg (for lobes com- 
pletely filled by relativistic plasma at the projected distance from 
the AGN, see Section [63] ) and the radiative timescale determined 
from the radio observations (170 Myr; GVM08) we can estimate 
the AGN power output to be 3.04x 10*^ erg s"'^. The required effi- 
ciency of the AGN in converting the co oled gas int o energy would 
thus be Econv = Prr,„rh / A/c^ =8 x 10~ ^ . ISuH ct alj (l2007) note that 
mass loss from stars may contribute to the mass of cool gas avail- 
able, reducing the required efficiency further. We estimate the stel- 
lar mass loss from AGB stars w ithin the core to be 0.03 MQyr"'^, 
based on the mass loss rates of lAthev et al.l ( l2002h . This suggests 
that the stellar mass loss rate may be at least capable of balancing 
cooling losses from the corona. 

If the AGN is fuelled directly from the hot gas of the corona 
(i.e., without the requirement that the gas first cool out of the X- 
ray regime), a first approximation of the energy avail able can be 
obtain ed by assuming the Bondi accretion rate (e.g., lAUen et al] 
l2006h . The Bondi accretion radius Ra ~ GAIbh /ca, where Ca 
is the adiabatic sound speed (~536 km s^^) and Mbh the black 
hole mass, which we can estimate from the relationship between 
Mbh and the central stellar velocity dispersion ( Gebhardt et al.l 
l200(i) . Taking the averaged velocity dispersion for NGC 6051 from 
LEDaB (7=343 km s"\ we find Mbh=9.1x 10** Mq, and there- 
fore a Bondi accretion radius i?A=13.65 pc. The Bondi accretion 
rate is defined as A4b ~ 'iTvR\p{RA)ca where p{Ra) is the den- 
sity at the Bondi accretion radius. Since this radius is a factor > 100 
smaller than our innermost spectral bin, it is clear that the uncer- 
tainty on p{Ra) will be the dominant source of uncertainty in our 
estimate of the accretion rate. 

We can place a lower limit on p{Ra) by assuming that it is 
equal to the mean density of the corona. Unlike cooling time and 
mass deposition rate, the Bondi accretion rate will therefore de- 
pend on the volume we assume the corona to occupy. Adopting 
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the lower density derived from the spectral extraction region we 
find Mb > 0.002 MQyr"^, implying a maximum conversion effi- 
ciency of Cconv < 0. 148, rather higher than the efficiency of 0. 1 of- 
ten assumed. If we instead estimate p{Ra) by fitting two /^-models 
to the deprojected density profile, we find a central electron density 
of Ue ~ 1.2 cm"'', about ten times greater than the mean den- 
sity of the central spectral bin. This leads to a Bondi accretion rate 
Mb = 0.022 Mqyi-^ and a conversion efficiency ficonu = 0.014. 
This is likely still an underestimate, as the 0.3-1 keV surface bright- 
ness profile, which most closely reflects the density of the cool gas, 
does not flatten in the core at the resolution of Chandra. However, 
these values are sufficient to confirm that accretion directly from 
the hot phase could fuel the AGN. We therefore conclude that at 
least in principle, the corona can support the AGN outburst over 
long periods without the need to build up significant quantities of 
cold gas. 

Given the strong temperature difference between the corona 
and ICM, conduction could have a significant effect on the corona 
gas. Studies of coronae in clusters across a range of temperatures 
suggest that conduction must be suppressed (by factors 30-50 0) for 
the corona to survive more than a few Myr dSun et alj|2007l) . This 
is an important factor in differentiating coronae from simple cool 
cores, as the suppression of conduction implies a physical separa- 
tion between corona and ICM, with the corona sustained by stellar 
mass loss within the galaxy rather than cooling from the cluster. 

We assume a temperature gradient of 1.36 keV across a dis- 
tance of 5.16 kpc, the distance between the midpoints of the first 
and second bins of our temperature profile. We find the mean free 
path of electrons in the IC M to be 170 p c. Conductivity follows 
the relations described by jSDitzei|[l962h where mean free path 
Ae <^ r/|dr/dr|, the scale height of the t emperature gradient, an d 
saturates where the two are comparable fcowie & McKee' ll977h . 
We find Ae ~ 185 pc and estimate the temperature gradient ex- 
tends over 5 kpc, so calcul ate the "classical" con duction rate (using 
the methods described in lO'SuUivan et alj2007l) . 

We find the classical rate of energy transfer to be 
^lO*^ erg s~^, assuming a radius of 1875 pc for the corona. 
The bolometric X-ray luminosity of the corona gas is only 
~5.7 X IC*" erg s^^, and radiative cooling therefore cannot balance 
conduction. The implied timescale for the corona to be heated to 
the temperature of the surrounding IGM is ~10 Myr. If we use the 
smaller radius estimated from the pressure balance argument, the 
conduction rate is reduced, but the corona will still be heated in 
<20 Myr. Taking the radiative age of ^--^170 Myr as the timescale 
over which the temperature gradient has existed, conduction into 
the corona must be suppressed by at least a factor 15-20, presum- 
ably owing to the effects of magnetic fields. Shorter timescales will 
reduce the need for suppression, but are unlikely to resolve the is- 
sue. 

The estimated mechanical power of the AGN jets would also 
heat the corona on short timescales if the jets were interacting with 
the corona gas. Assuming a jet power of ~ 10*"^ erg s^^ (see Sec- 
tion |6.5l ) only 0.4 per cent efficiency would be required to heat the 
corona over the timescale of the current outburst (170 Myr). The 
apparent tight collimation of the jets within the corona suggests that 
any interaction with the corona gas is minimal. We conclude that 
neither conduction nor the AGN jets are heating the corona to a sig- 
nificant degree. This suggests that the presence of the corona may 
disrupt the feedback relationship between gas cooling and AGN 
heating, since the jets do not reheat the material which fuels the 
nuclear activity. 
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Figure 7. GMRT 610 MHz image of 4C+24.36, with regions used to de- 
rive radio pressures (selected based on the 235 MHz image) marked by 
ellipses. The white cross marks the position of the radio core, at the centre 
ofNGC6051. 

6 PROPERTIES OF THE RADIO SOURCE AND IMPACT 
ON THE ICM 

6.1 Pressure balance 

The absence of strong X-ray surface brightness features associated 
with the radio source suggests a) that the jets and lobes of 4C+24.36 
are in approximate pressure equilibrium with the ICM, and b) that 
any cavities associated with the lobes are only partially filled by 
radio plasma, or that some other factor is affecting our estimate of 
the expected surface brightness deficit. A comparison of the ap- 
parent pressures of the thermal and relativistic pla sma can provide 
insight into the particle content of the jets (e.g., iDunn & FabianI 
l2004l ; lDunn et alJl2005l : lBirzan et al . 2008). We therefore estimate 
pressure profiles from both the X-ray and radio data. 

Our radio analysis is performed using the GMRT 235, 327 
and 610 MHz observations presented in GVM08 and archival 1.4 
and 4.9 GHz VLA data. We select four regions along each jet, with 
minimum size based on the lowest resolution data (~13" HPBW at 
235 MHz). These regions are shown in Figurel?] While they appear 
large compared to the structures in the 610 MHz image, they are 
well matched to the size of jets and lobes at 235 MHz and therefore 
suitable for examination of the spectral index across the different 
frequencies. The volumes of the regions are estimated assuming 
them to be oblate or prolate ellipsoids with rotation axes aligned 
with the jet. Fluxes are measured in each band to obtain the inte- 
grated radio spectra for each regio n. The spectra are fitted with the 
Synage-l~l- package ('MurHia"200l'), adopting a continuous injection 
model (CI, Kardashev 1962) to derive the injection spectral index 
Olinj for each region. Using Qi„j , we then derive the physical pa- 
rameters in each regions assuming minimum energy conditions, in 
which the contribution to the total energy content from relativistic 
particles and magnetic field are approximately equal. Uncertainties 
on derived parameters are estimated using the errors on (y.inj pro- 
vided by the best-fits, and an assumed 5% uncertainty on the radio 
flux. For our analysis, the synchrotron radio pressure is defined to 
be 

,Up BL„ {l + k)E, 

Psync = f/fl + — = 3^_^ , (3) 

where Ub and Up are the energy density of the magnetic field 
and relativistic particles, Bmin and Ee are the minimum energy 
magnetic field and total energy in electrons, is the permeability 



of free space (47r x 10~^ NA~^ or Atv if working in Gauss and 
cgs units), V is the volume of the region, (j) is the filling factor, and 
k is the ratio of energy in non-radiating particles to the energy in 
electrons. 

For the initial pressure calculations we assume (j) — 1 and 
k = 1, implying that half of the energy in particles is in the form of 
non-radiating particles, as would be the case in an electron-proton 
jet. An electron-positron jet would have k = 0. We note that our 
definition of k dif fers from some commonly used formulae (e.g., 
iFabian et alj|2002h which use the ratio of the total energy in parti- 
cles to the energy in electrons (i.e. Ep = kpE^). Under our defi- 
nition, = — 1, so our assumption of equal energy in radiating 
and non-radiating particles would imply kp = 2. 

The pressure of the relativistic plasma also depends strongly 
on the range of energies in the particle population. As particles with 
low 7 values radiate at frequencies which are not practically ob- 
servable in most sources, it is generally necessary to assume an 
energy range. We have chosen to consider three possible ranges: 

(i) 7 values equivalent to the frequency range 10 MHz- 100 GHz, 
which we refer to as the standard equipartition case. This is com- 
parable to many previous studies in the literature. 

(ii) 7=100-5500, referred to as revised equipartition with 
'y^i„=100. This is the range chosen by GVM08. 

(iii) 7=10-5500, referred to as revised equipartition with 
7„i„=iO, chosen to allow us to examine the effect of including 
low energy electrons, under the assumption that the electron energy 
distribution continues to follow a power law down to low Lorentz 
factors. 

We discuss these ranges and the reason for their choice below. 

Most previous studies of the physical parameters of the jets 
and lobes of radio galaxies have adopted the frequency range of 10 
MHz- 100 GHz normall y used in the standard equipartition equa- 
tions jPacholczvk " 1970). In many cases high k values (of the order 
of 100-10000, e.g., Dunn et al. 2005) are found, indicating a re- 
quirement for a large fraction of energy in non-radiating particles. 
Inverse Compton observations of high-redshift radio sources have 
shown that the assumption of a power law distribution of electron 
energies is reliable down to 7 values of a few hundred, which would 
have synchrotron frequencies ~ 500 kHz. These lower energy elec- 
trons could provide at least part o f the additional pres s ure implied 
by th e large observed k values l lCroston et al.l 1200 Sl : iDunn et afl 
I2OIOI) . GVM08 include these lower energy electrons in their es- 
timates of total energy an d pressure, using the revised equiparti- 
tion equations of Brunetti et al.l j 19970 with a minimum electron 
energy cutoff (7™,™) instead of a minimum frequency. Following 
GVM08, we estimate the radio pressure assuming 7min ~ 100 and 
a ~ ctinj- The maximum energy cutoff, 'ymax, was selected such 
that the maximum emitting frequency of the electrons matched the 
observed break frequency of the radio spectrum, again following 
GVM08. In the lobes, this gave 7maa;=5500. Approximate break 
frequencies for the different regions are given in Table[3] 

The thermal pressure of the ICM was derived from the depro- 
jected spectral profiles with annuli selected to match those of the 
ellipses used in the radio. The pressure in each annulus was calcu- 
lated as Pth = nkpT where we have assumed an ideal gas with 
n = 2ne . The resulting profile agrees well with the pressure pro- 
file derived from the XMM-Newton analysis of OS05, except in the 
central bin which the XMM profile does not resolve. We note that 
the X-ray pressure profile presented in GVM08 was unfortunately 
incorrect, with the normalisation of the thermal gas pressure profile 
decreased by a factor of 10 from its true value. Figure [8] shows a 
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Figure 8. Radial profiles of X-ray derived thermal gas pressure and equipartition synchrotron pressure for the jets and lobes of 4C+24.36. Black crosses 
indicate the thermal pressure derived from deprojected Chandra spectral models, the black line and grey shaded region the XMM pressure profile and Icr 
unceilainty. The red crosses indicate the synchrotron pressure under our assumption of the revised equipartition conditions with 7min=100, and the green 
diamonds the pressure in the lobes if 7mm=10. The blue stars represent the "standard equipartition" pressure estimates for the lobes. Vertical symbol sizes 
indicate la uncertainties for the thermal and revised equipartition synchrotron pressures. 



Table 3. Break frequencies measured in each region of the jets and lobes 



Region 


f break 


Region 


f break 




(GHz) 




(GHz) 


El 


30 


Wl 


5 


E2 


3 


Wl 


2 


E3 


1.2 


W3 


1 


E4 


0.7 


W4 


0.7 



Note that in regions Wl and El the frequencies are poorly constrained, as 
the change in powerlaw slope is small. 

comparison between the thermal gas pressure profiles and the min- 
imum synchrotron pressure estimates described above. 

The estimated synchrotron pressure in the jets and lobes is 
lower than the thermal pressure at all radii. This is expected; pres- 
sure imbalances between radio sources and the surrounding ICM 
are common in extended FR-I radio sources, with the ratio of 
thermal to synchrotron pressures (PthfPsync) having values of up 
to ~ 100 (under standard equipartition assumptions, iFeretti et al.l 
Il992l) . The uncertainties in the synchrotron pressures in the outer 
bins are large. For our "/min ~ 100 estimates, the imbalance is 
3 — 5(T significant, except in the east lobe where the significance 
drops to 2.2(7. While this means that statistically the east lobe may 
be in pressure equilibrium, the consistently low pressure estimates 
in the other regions strongly indicate an imbalance in the source as 
a whole. 

The synchrotron pressures decline by a factor 3-4 with radius, 



as expected, given radiative energy losses from the particle popu- 
lation of the jets. The pressure profiles are roughly consistent with 
a constant pressure ratio of Pth/Psync ~15. There is no indication 
of a strong gradient in pressure ratio, which might have been ex- 
pected if the jets are entraining significant quantities of ICM gas 
over their whole length. 

As our assumed values of k=l and (j>=l appear to be incon- 
sistent with the data for the 7mi„ — 100 case, we estimate the 
change in these parameters which would be necessary to produce 
pressure equilibrium in the lobes. Setting equation[3]to be equal to 
the thermal pressure and rearranging allows us to place limits on 
these parameters, shown in Table |4] 

Our alternative estimates of pressure in the lobes can also be 
considered in this way. Assuming standard equipartition conditions 
decreases the estimated pressure in the lobes by a large factor, with 
pressure ratios of ~ 300 — 350. This translates into a high value of 
( 1 + fc) /(/), which strongly suggests a large amount of additional en- 
ergy in the non-radiating particle population. The revised equipar- 
tition estimates for "/^i„=10 reduce the pressure imbalance to a 
factor of 3.9-4.5, with Icr uncertainties consistent with the thermal 
pressure. While these uncertainties are large, owing to the extrap- 
olation of uncertainties in the spectral index to low 7 values, this 
indicates that the additional energy required for pressure balance 
could be provided by the inclusion of electrons with low Lorentz 
factors. The remaining imbalance could be explained by relatively 
minor changes in k or the latter is of particular interest given 
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Table 4. Magnetic field, pressure, filling factor and particle content of the radio lobes (regions E4 and W4) 



Lobe 


Pth 
(erg cm~3) 


revised equipartition, 7min 

Bmin Psync 

ipG) (erg cm- 3) 


= 100 

(l + fc)/«i 




standard equipartition 

Psync (1 + k)/(j, 

(erg cm- 3) 


revised equipartition, 7min = 10 

Bmin Psync {1 + k) / <j) 

(p.G) (erg cm-3) 


East 
West 


31.40±0.98 X 10'" 
24.80±1.00 X 10-" 


6.36+^-?^ 2M+fll X 10-12 

emtlil i-96l?;?S X 10-" 


160.3l-3.93 

123 n+135.2 


1.92 
1.85 


9.11x10-1* 3192.6 
8.44x10-1-1 2448.2 


11.45 6.95x10-12 
10.97 6.38x10-12 


43.3 
32.4 



Note that values of (1 + fc)/0 are those required for pressure equilibrium in the three cases. The spectral index used was di^j = 1.01. 



the weak cavities in AWM 4. The large differences between these 
two estimates clearly indicates the importance of the equipartition 
assumptions in determining the pressure imbalance and any inter- 
pretation of the data. 



6.2 Outburst Timescale 

The activity timescal e of the radio sourc e can be estimated in sev- 
eral ways. FoUowing lParma et alj jl983) . GVM08 estimate the age 
of the source by modeUing the change in radio spectral index along 
the jets and lobes, under the assumption that the age of the ra- 
dio emitting electrons increases linearly with distance. They find 
ages of 17ll^g Myr and I59I22 ^^^^ west jets and 

lobes respectively, assuming 7mm=10 and a fitted injection spec- 
tral index ai„j ~0.5. Recalculating with 7mi„=100, we find ages 
of 183124 Myr and I64I23 Myr for the east and west jets and lobes. 
As these calculations are based on multiple independent spectral in- 
dex measurements in each direction, this is likely the most reliable 
radiative age measurement available. However, it should be noted 
that the ages are most applicable to the lobes; a younger age would 
be found at a given point along the jet. GVM08 also note that the 
estimated age is dependent on the assumption that losses associ- 
ated with source expansion are negligible. If this is not that case, 
the radiative age could be reduced by up to a factor of ~3. 

Radiative ages can also be estimated for individual regions, 
based on the minimum energy magnetic field measurements, us- 
ing Equation 1 of GVM08. The injection spectral index estimated 
for the lobes (a ~ 1) differs from the global estimate (Qmj ~ 
0.5), which can lead to large differences in radiative age. Adopt- 
ing 7mi„=100, we find an age of 1351^^ Myr (1321^2 j^yj.) fo^. 
the east (west) lobe. These values are consistent with those de- 
rived from the ages derived from modelling of spectral aging. 
The standard equipartition assumptions produce much greater ages, 
~252 Myr and ~236 Myr for the east and west lobes respectively. 
Conversely, adopting the revised equipartition assumptions with 
'ymin=W, we find much shorter ages, --^66 Myr for both lobes. 

Independent limits on the age of the radio source, and there- 
fore on its magnetic field strength and particle content, can be esti- 
mated based on dynamical arguments. A lower limit on the age of 
the lobes can be estimated from the time taken for them to grow to 
their current size assuming that (as we see no evidence of shocks 
associated with them) their expansion has been subsonic. We can 
also estimate the time taken for the jets to expand to their ob- 
served length at the sound speed (~770 km s^^ for fcsr=2.6 keV), 
though there is a possibility that the jets may have expanded super- 
sonically earlier in their history and the shocks produced have now 
moved out of the field of view or weakened to become undetectable. 

A more realistic value may be the time taken for the lobes to 
rise buoyantly to their current position. The buoyant velocity is 



where R is the mean radius from the core at which the lobe 
is found, S is its cross-sectional area in the direction of motion, V 
its volume, AI{< R) the total gravitational mass w ithin R, and Cd 
the co efficient of drag, typically taken to be 0.75 jChurazov et al.l 
I2OOII) . Finally, we can also estimate the time required for the 
ICM to refill the displaced volume as the lobe rises, trefui = 
2Ry'r/GM{< R), where r is the mean radial size of the lobe. 
These four timescale estimates are given in table[5] 

Using these timescales to estimate the magnetic field strength, 
we can again assume pressure equilibrium so as to determine 
(1 + fc) /(j}. However, we must take into account the dependence of 
the energy of Jheparticle p opula tion on the magnetic field strength. 
Following ISrunetti et al.l l ll997l) , this dependence is of the form 
Ep oc B^"*"", leading us to modify equation[3]as 



1 + k 



Pth 



2fio 



W 
C 



B 



(l + a) 



(5) 



Vbu 



/2GM(< R)V 



SR^Cd 



(4) 



where C is a constant related to the synchrotron flux, -jmin 
and the spectral index. We can therefore calculate the change in 
{1 + k)/(j) relative to our original value. The timescales and result- 
ing {1 + k)/(l) values, assuming magnetic field strengths derived 
for the revised equipartition conditions and ^^i„=\00, are given in 
Table [5] Estimates based on the sonic expansion timescale of the 
lobes are not included as this is short enough to imply synchrotron 
pressures greater than the thermal pressure. 

The dynamical timescales are all shorter than the GVM08 ra- 
diative age derived from the spectral aging along the jets. The buoy- 
ancy and refill timescales of the west lobe are comparable to the 
radiative age derived from the 7„i„=100 magnetic field measure- 
ment in that lobe. However, both timescales depend strongly on 
the filling factor; reducing would produce a shorter refill time, 
but a longer buoyant rise time. The change in buoyancy timescale 
will depend on the morphology of the lobe material. For a simple 
spherical cavity, ttuoy for (j)=0.2, tbuoy will increase by 

a factor 1 .3 to ~ 175 Myr for the west lobe. The clumpy appearance 
of the radio lobes suggests that they have a large surface area for 
their volume, so tbuoy is likely to increase more slowly. However, 
it seems possible that filling factor considerations account for part 
of the difference between the buoyancy and radiative timescales. 

The sonic expansion timescale of the jets is comparable to the 
'y„j„=10 radiative age in the east jet, and longer than it in the west 
jet. The large uncertainties on the radiative age mean that there is 
no formal disagreement, but this suggests that the jets expanded 
supersonically over some portion of their history. If this is the case, 
shocks would have been produced, perhaps providing additional 
heating. 



6.3 Filling factors and Mixing in the Lobes 

The lack of clearly detected cavities in AWM 4 raises the possibility 
that the lobes are not completely filled by relativistic plasma. Some 
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Table 5. Radiative and dynamical timescale estimates. 



Lobe 


trad, 
Imin = 100 

(Myr) 
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(Myr) 


trad 
Imin = 100 
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^min — 10 
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ts,lohe 

(Myr) 


^s,jet 
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^buoy 

(Myr) 


trefUl (1 - 

(Myr) 


f k)/(j)s,jst (1 - 


t- k)/<t>buov 


(1 + k)/(f)refM 


East 
West 




171+« 


135+^*' 

1 17+72 
"^-52 


/t£+47 
66I41 


27.53 
26.77 


69.5 
87.3 


103.16 
134.76 


120.02 
129.20 


430.3 
239.0 


249.6 
118.8 


196.3 
128.7 



trad,total is the radiative age estimated from spectral aging along the jets 

trad lobe is the radiative age estimated from the spectral index of the lobes 

ts lobe is the time required for lobes to expand to their current width at the sound speed 

tsjet is the time required for the west jet and lobe to reach their current length at the sound speed. 



support for this provided by the radio maps of 4C+24.36 (e.g., Fig- 
ure |7), which show both lobes to be clumpy, with brighter regions 
apparently tracing the jets out to the tips of the lobes and fainter, 
flocculent regions to north and south. By comparison, the jets are 
much more sharply defined. 

The lobes may have formed because the jets become uncol- 
limated, with the relativistic plasma breaking up to form smaller 
clouds and filaments which then mix with the ICM. The appear- 
ance of the radio emission and the lack of evidence for higher ICM 
temperatures along the lines of sight to the lobes argue against di- 
rect mixing of the jet plasma with the ICM gas, instead suggesting 
that it remains confined, probably in the form of clumps with em- 
bedded magnetic fields. For this scenario, we can place limits on 
the filling factor based on the X-ray surface brightness in the 
lobes. 

For the western lobe, where only a weak deficit in surface 
brightness is seen, we estimate a 3a upper limit of (f} < 0.76, as- 
suming that there are no significant additional sources of X-ray flux 
in the lobe. The observed decrement of ~40 counts (1120 counts 
expected, 1080 observed, 1-3 keV) suggests 0=0.21, and is consis- 
tent with a zero filling factor. As discussed in Section[3] we would 
expect a deficit of 184 counts if the lobe were empty of thermal 
plasma. 

In the east lobe, where we see some indication of a cavity, 
we use the circular region shown in Figure [2] (as compared to the 
larger ellipse discussed in Section |3]l and find ~430 counts, giv- 
ing a decrement of ~60 counts below the ~495 expected from the 
surface brightness model. Assuming ICM properties typical at this 
radius, this implies a that only ^^24% of the volume of the lobe is 
filled by relativistic plasma. However, the X-ray emission over the 
remainder of the lobe is uneven, with weak filamentary structures 
to north and south of the holes, so our estimate of the filling fac- 
tor may be too low. Considering the lobe region as a whole (i.e., 
the ellipse region in Figure |2]l, we can place a 3a upper limit of 
(I) < 0.43. 

As discussed in Section [631 Psync oc provided the mag- 
netic field B is similar to the equipartition or minimum energy field 
strengths. Using the revised equipartition conditions, the magnetic 
field can be defined as 

„ f 67r(l + fc)ei2Ls^nc V^'^^"' 

^-=(, ) ' 

where Lsy„c is the sy nchrotron luminosity, V is the volume, 
and a is the spectral index jGovoni & Ferettill2004l) . The constant 
C12 is a function of the spectral index and the frequency band over 
which we integrate ( |Pacholczvklll970l) . 

For our source a ~ 1, so Payne oc ^ (1 + k) / (p. Our filling 
factors thus suggest a pressure increase of a factor of 2-2.2, with 
a 3a minimum increase of ^^1.5. This is insufficient to resolve the 



pressure imbalance we estimate from the revised equipartition as- 
sumptions with 7mi„=100, but would significantly reduce the dis- 
crepancy. Based on these limits on <f) we can estimate the value of 
k (and the 3a upper limits) to be k=31.5 (k < 741.6) for the east 
lobe, and A:=24.8 (k < 517.0) for the west lobe. Filling factors 
of order <^=0.015 would be required to match the range of pres- 
sure imbalance in the lobes. While this would be consistent with 
the lack of surface brightness features in the west lobe, it is prob- 
lematic for the decrement seen in the east lobe. This suggests that 
simple mixing with the ICM gas cannot account for the observed 
pressure imbalance. 

If we assume 7mi„=10, the limits on k will of course be re- 
duced. In this case A;=10.4 for the east and fc=6.8 for the west lobe. 
Within the uncertainties, the "f^i„=lO pressures are consistent with 
the k=l expected for an electron/proton jet with equal energies in 
each class of particles. However, the uncertainties on k are very 
large, and would also be consistent with fc=0 or with the values 
estimated using 7„i„=100. 

6.4 Entrainment in the Jets 

A possible additional source of pressure is gas entrained and heated 
by the jets. It is thought tha t entrainment slows FR-I jets and 
causes their broadening (e.g.. IWorrallll200^ , and heating of en- 
trained material is likely to be most effective where the jet is still 
coUimated and its velocity is high. As we see little sign of inter- 
action between the jets and corona, the most likely source of en- 
trained material is mass loss from stars inside the kpc-scale jets. 
We place a maximum limit on the jet width from the 0.4" res- 
olution of the 4.9 GHz data (250 pc), and take the length from 
mean distance between the radio core and the points where the 
jets begin to flare (1.5 kpc). We estimate the stellar mass loss rate 
for stars within the jets from the mass loss rate of AGB stars in 
old stellar populat ions, M. = 0.078(Lb/10"Lso) Mq yr~^ 
jAthev et alj|2002l) . using the mean stellar density within 1.5 kpc, 
2.76xlO*I/B0 kpc~^. Each jet is therefore likely to entrain ap- 
proximately 1.63x10""* A/0 yr"\ If we assume the jet has been 
active for 170 Myr and is rapid enough to transfer all entrained ma- 
terial into the lobes, we expect ~ 2.8 x 10* Mq of stellar material 
to be mixed with the relativistic plasma. 

Based on the ICM density profile and the ellipsoidal lobe re- 
gions assumed above, we estimate that if they are completely filled 
by radio plasma, the lobes have each displaced ~ 3.7 x 10^ Mq of 
ICM gas. The pressure imbalance suggests that the radio compo- 
nent in fact occupies only <25 per cent of the lobe volume. For en- 
trained stellar wind gas to provide the remaining pressure, it would 
need to be heated to temperatures >5 MeV. Such material could be 
the source of additional particles in the jet implied by fc > 1, with- 
out contributing significantly to the X-ray flux from the lobes. It 
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Figure 9. The variation in instantaneous loss timescale with density, with 
each solid diagonal Hne representing particles of a given 7 factor as labelled 
on the plot. A characteristic age of 170 Myr is shown for reference. 



would produce inverse-Compton scattering, but an electron energy 
of 5 MeV is equivalent to 7e ~ 10, so CMB photons will be scat- 
tered into the mid-Infrared rather than the X-ray band. Assuming 
that this population provides all of the pressure required to bring the 
lobes into equilibrium, and a nominal spread of energies of a factor 
of four, we can expect an Infrared flux of ^1.1x10"^" erg s~^. 
This is well below most estim ates of the cosmic Infrared back- 
ground (e.g.. lDoleet"ai]|2006l) and approximate ly three orders of 
magnitude below the 15/xm limit from ISOCAM telbazetaljl"999l) 
which most closely matches the expected energy. 

It is also possible that ICM and stellar material are entrained 
and heated in the broader section of the jets which extends from 
the corona out to the lobes. Since the velocity of the jet is expected 
to be lower, and the broader jet is likely to entrain more material, 
it is possible that the heating could be mild, producing a separate 
high temperature thermal component rather than adding to the jet 
material. However, if the lobes are in equilibrium, the maximum 
energy involved in heating entrained gas is that required to achieve 
pressure equilibrium in the volume filled by jet material. Assum- 
ing (/!>=0.25, this is ~ 5 X 10^® erg, a significant fraction of the 
mechanical energy of the jet. 

A final consideration is the effect of material entrained within 
the jet on the relativistic particle population. Coulomb losses due 
to collisions with thermal electrons could be important for particles 
with l ow 7 factors, if the two plasmas are not segregated. ISarazii] 
( Il999h estimates the loss rate to be 



1.2 X lO^^^n, 



1.0 + 



ln{')/ne 



75 



(7) 



The instantaneous loss timescale tioas ~ 1 /bcouiil) then 
gives an approximate timescale over which particles will lose their 
energy via Coulomb interactions. Figure [9] shows the relationship 
between loss timescale and density for particles of different 7 fac- 
tor. 

From the estimated mass of gas lost from stars in the kpc- 
scale jets, we can estimate the maximum density of thermal plasma 



within the jets in this region. Taking a lower limit on the jet ve- 
locity from the large-scale expansion speed of the radio source 
(0.002c, GVM08), we estimate a maximum density of rig ~ 3 x 
10~* cm"''. This suggests that Coulomb losses could be signifi- 
cant for particles with 7 < 10 but that more energetic particles will 
not lose a significant amount of energy. If stellar mass loss is the 
only source of entrained material, losses in the lobes will be even 
smaller, as the greater mass of gas occupies a much larger volume. 

However, it is clear that if the entrained material is signifi- 
cantly denser. Coulomb losses will rapidly become important. If 
densities > 10"'^ cm~^occur along much of the jet, we would ex- 
pect significant energy losses from particles with 7 < 100. In this 
case the entrained and heated gas would have to provide the nec- 
essary additional pressure, since the contribution from low-7 rela- 
tivistic particles would be removed. 

We can rule out the most extreme case, entrainment of un- 
healed (or only mildly heated) ICM plasma by the jets. For the ob- 
served ICM densities. Coulomb losses would affect particles radi- 
ating at radio frequencies, and would produce a clear change in the 
radio spectrum. This could only be avoided if the thermal plasma 
were segregated from the relativistic particles, with the thermal 
plasma entrained in non-interacting clumps. 

In summary, the lack of clearly detected cavities and the 
clumpy, filamentary appearance of the radio lobes suggest that the 
lobes are regions of mixing between the relativistic jet plasma and 
the ICM. Our best estimate of the filling factor is (/)=0.2-0.25. This 
naturally leads to a reduction in estimates of the energy in non- 
radiating particles. If the energy distribution of the electron pop- 
ulation follows the measured powerlaw to 7 = 10, then no ad- 
ditional components are required to bring the radio plasma into 
pressure equilibrium with its surroundings, albeit with large uncer- 
tainties. Alternatively, additional material could be entrained by the 
jets, probably without producing detectable extra emission compo- 
nents. However, heating such material to the temperatures required 
to produce pressure equilibrium would require a significant frac- 
tion of the total energy of the jets. Coulomb interactions between 
entrained material and the relativistic plasma of the jet could also be 
important, and while material from stars within the kpc-scale jets 
is unlikely to significantly affect the particle population, entrain- 
ment of larger quantities of gas from the ICM could cause the loss 
of the lowest 7 electrons, and any contribution toward the pressure 
support of the lobes they provide. 



6.5 Energetics of the radio source 

To obtain a first approximation of the maximum mechanical en- 
ergy output of 4C-I- 24.36, we can assume that the radio lobes have 
in fact excavated cavities in the ICM, and estimate the power re- 
quired. The enthalpy in a cavity is given by _ff = 'yPV/('y — 1), 
where 7 is the ratio of specific heats (5/3 for a non-relativistic gas, 
4/3 for relativistic gas). Assuming the cavities to be completely 
filled with relativistic synchrotron-emitting plasma, the mechan- 
ical energy required is 4PV. Assuming the projected distances 
of the lobes and using the ellipsoidal regions described in Sec- 
tion [3] (and shown in outline in Figure O, we estimate this to be 
1.63±0.02 X 10^" erg. Taking the age estimate of 170 Myr from 
radio spectral modelling, the power output of the AGN would be 
3.04±0.04 X 10*^ ergs"^. The bolometric X-ray luminosity of 
the ICM within a sphere with radius equal to the outer radius of the 
western lobe is ~ 1.29 x 10*"^ erg s~^, neglecting the contribution 
from the corona. 

At first glance, the similarity between the energy available 
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from the AGN and the radiative energy loss suggests a balance 
between heating and cooling. OS05 discussed the possibility that 
AWM 4 may at one time have possessed a large cool core com- 
parable to those observed in clusters of similar temperature and 
galaxy population. The estimated that heating such a cool core to 
the present typical temperature of 2.6 keV would require roughly 
9 X 10^* erg. At best only about half of the AGN energy input would 
be available to heat the gas after radiative losses, but over the life- 
time of the outburst, the total energy available is ~ 9.4 x 10^* erg, 
a close match to energy required. 

This simple calculation may underestimate the energy sup- 
plied by the AGN. The timescale of activity could be shorter, im- 
plying a greater power output for the same total energy and smaller 
fractional losses from radiative cooling. For a timescale of 1 35 Myr, 
similar to the buoyant rise time of the west lobe, the power of the 
jets would be ~ 3.8 x lO"^ erg s"\ 

As discussed above, there is also the possibility that the jets 
may have driven shocks into the surrounding ICM during earlier 
phases of the outburst. Despite the difficulty of observing such 
shocks, owing to the short period between their formation and 
their progression to the low density outer regions of the cluster 
halo, several examples have been confirm e d (e.g..[Gitt i et al. 201Q; 
Graham et all |2008| ; iForman et all l2007l : iMcNamara et alj bood : 
Nulsen et alj2005l) . These tend to be weak shocks with Mach num- 
bers M < 1.4, but may transfer considerable amounts of energy 
into the ICM. The shocks in the Perseus and Virgo clusters involve 
considerably more energy than the radio lobe cavities in those sys- 
tems, sug gesting that shocks a re the dominant mechanism of energy 
injection jporman et alj|200^ . and simulations of jet/lCM interac- 
tions show that multi ple shocks may be produced during a single 
nuclear activity cycle (Briiggen et al. 2007). It is therefore plausible 
that shocks have contributed to the heating of the ICM in AWM 4 
at some stage, though none are currently detected. 

However, it seems more likely that our calculation overesti- 
mates the energy input from the AGN. As discussed in Section [63] 
the lobes are probably partially filled by relativistic plasma, im- 
plying that we have overestimated the cavity volume. To examine 
the effect of the filling factors calculated in section [63] we con- 
sider the lower estimate of filling factor (j> — 0.21 and the upper 
limit of (j) < 0.43. These imply enthalpies of 3.4x10^* erg (< 
7.0 X 10®* erg) and mechanical jet powers of 0.64x10''^ erg s"^ 
(< 1.31 X 10^'^ erg s~^), assuming a simple scaling with total vol- 
ume. This suggests that the energy available from the radio source 
is just sufficient to balance cooling, unless some other mechanism 
for energy injection (e.g., shock heating) plays a significant role. 



6.6 Non-thermal X-ray emission from the radio lobes 

An alternative reason for the lack of clear cavities could be the pres- 
ence of non-thermal X-ray emission from the radio lobes, which 
would "fill in" some or all of the cavity surface brightness decre- 
ment. This motivates us to determine the limits of any inverse- 
Compton component. 

Of the two sources of seed photons available, the energy den- 
sity of radio-frequency synchrotron photons produced in the lobes 
is a factor ^ 10~® lower than that of cosmic microwave back- 
ground (CMB) photons, and we therefore neglect synchrotron self- 
Compton scattering. We can estimate the expected X-ray flux from 
inverse-Compton scattering of CMB photons from the radio spec- 
tral properties, under the assumption of equipartition. We can de- 
termine which electrons are involved in scattering from the relation 
between the change in photon energy and the Lorentz factor of the 



Table 6. Expected and measured non-thermal X-ray emission in the lobes 







East lobe 


West lobe 


a 




1.01 


1.01 


Eparticles (10^'^ erg) 




2.088 


2.021 


expected Fx ( erg s~^. 


0.7-7 keV) 


1.38xl0~i5 


1.33x10-1^ 


measured Fx ( erg s~^ 


, 0.7-7 keV) 


<8.21xl0-i5 


<4.62x 10-15 



electrons, irx/vcMS — (4/3)7^, where i^x and vqmb are the fre- 
quencies of the scattered X-ray photon and the pre-scattering CMB 
photon, and 7e is the Lorentz factor of the scattering electron. To 
scatter a photon at the peak frequency of the CMB, ~160 GHz, 
up to an energy of 1 keV, an electron with 7^ ~ 1000 is required 
(7e ~ 330 for 0.1 keV, ~2800 for 7 keV). This is within the range 
of L orentz factors used i n our equipartition calculations. 

lErlund et al. I l l2006h provide an estimate of the approximate 
inverse-Compton (IC) luminosity in a given energy band, under the 
simplifying assumption of a monochromatic distribution of CMB 
photons: 



Lx = — - — Nootc — T— 



1 



(8) 



where a is the spectral index determined from the radio data, 
(TT is the Thompson scattering cross-section, c is the speed of light, 
o-TcMB^Ucm^, the energy density of the cosmic microwave back- 
ground, and v\ and V2 are the lower and upper frequency bounds 
of the chosen energy range. The energy distribution of the electron 
population is assumed to be a power law, described by: 



E ■- 



iVo7'" 



^nieC d7, 



(9) 



where E is the energy of the electron population between 
Lorentz factors 'ymax and 7mm, A'^o is the normalisation of the 
power-law, is the electron mass, and p is the index of the power 
law, related to the spectral index of the radio emission and photon 
index of the X-ray emission by p = 2q -I- 1 = 2r — 1. 

Under the assumption of equipartition, the energy of the elec- 
tron population in the lobes can be determined from modelling of 
the radio spectral distribution. Table [6] lists the relevant radio pa- 
rameters and resulting flux estimates. We use the injection spectral 
index, as this is probably a more accurate estimate of the index at 
low values of 7^. We note that the energy in the relativistic particle 
population, Eparticles, IS Calculated assuming that electrons and 
positrons each contribute 50 per cent of the energy. The expected 
Fx is therefore an upper limit; for an electron/proton plasma the 
predicted flux will decrease by 50 per cent, and if a larger frac- 
tion of the jet energy is in heavy particles the flux will be further 
decreased. 

In order to place limits on the actual IC flux in the lobes, we 
extract X-ray spectra from each lobe and three annular bins ex- 
tending from the lobes out to our limiting radius of 395", and use 
a deprojection to model the contribution of gas along the line of 
sight to the lobes. The thermal emission is again modelled with 
an APEC plasma model, and we add a powerlaw model with fixed 
photon index 7 = a + 1 to account for the IC emission. Com- 
paring fits with and without the predicted levels of IC emission we 
find no significant differences in spectral parameters. The expected 
fluxes are a factor I0~^-I0~* of the measured total flux, and the 
photon indices of ~2 produce a similar spectral shape to the ob- 
served A:_Br~2.6 keV plasma spectrum. Allowing the normalisa- 
tion of the power-law component to fit freely, we find that the flux 
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is consistent with the predicted values and with zero, within 90 per 
cent uncertainties. The upper bounds on the IC flux are shown in 
table|6] Even assuming these maximum IC fluxes to be correct, the 
parameters of the plasma model are largely unaffected. We there- 
fore conclude that non-thermal emission is not detected from the 
lobes of 4C+24.36, and cannot be responsible for the lack of dis- 
tinct X-ray cavities associated with the lobes. 

6.7 Bending of the jets along the line of sight 

A final factor which could reduce the expected surface brightness 
decrement of cavities is bending of the jets. Bending would place 
the lobes at greater radii, where the density of the ICM is lower, 
reducing the pressure difference. While GVM08 note that both the 
small scale jet/counter-jet brightness ratio and the symmetry of jet 
properties on large scales suggest the radio source is aligned close 
to the plane of the sky, this does not preclude some bending of the 
jets along the line of sight. 

If we assume that the lobes rise buoyantly and the radiative 
age estimate of 170 Myr is correct, we expect the lobe will rise to a 
distance at which M{< R)/R'^ = 1.73 x 10^ Mq kpc""*. Using 
the gravitational mass profile of OS05, we find that this occurs at 
R ^ 77 kpc, compared to projected distances of ~53 kpc for the 
east lobe and ^^67 kpc for the west lobe. The thermal pressure at 
this radius is ~ 2 x 10^^^ erg cm"'', still a factor '--^lO greater 
than the synchrotron pressure in the lobes. This would suggest that 
bending could contribute ~20% of the apparent imbalance in the 
west lobe and ~35 per cent of the imbalance in the east lobe. Larger 
contributions imply higher velocities and pressure balance cannot 
be achieved within the radiative timescale as this would require 
spersonic motion. 

If the jets are bent in the line of sight, NGC 6051 must be 
in motion relative to the surrounding ICM. A true radius for the 
lobes of 77 kpc suggests a galaxy velocity of 3 16 km s"^ of which 
218 km s~^ would be along the line of sight. However, optical mea- 
surements show NGC 605 1 to have a velocity identical (within un- 
certainties) to the mean velocity of the cluster galaxy population 
jKoranvi & Gellej2002h . The uncertainties on these measurements 
are small enough that a velocity offset of the magnitude required 
appears unlikely. The large difference in luminosity (and hence 
mass) between NGC 605 1 and the other cluster galaxies, and lack 
of apparent substructure in the cluster also argues against recent 
interactions which could have produced such a velocity offset. We 
therefore conclude that NGC 605 1 probably has only a small veloc- 
ity along the line of sight, relative to the cluster, and that bending 
makes only a small contribution to the apparent pressure imbalance. 



7 DISCUSSION 

The analysis described above provides solutions to a number of the 
outstanding questions raised by previous observations of AWM 4, 
and helps provide a more coherent context for the cluster among 
other systems with powerful central AGN. However, a number of 
issues remain unresolved. 

While there are a number of weak features suggesting that the 
radio source 4C-I- 24.36 is interacting with the ICM, few are statisti- 
cally significant even in this relatively deep Chandra pointing. The 
lack of strong cavities associated with the radio lobes suggests that 
relativistic plasma only partially fills these volumes, mixing with 
ICM gas which may occupy as much as 80 per cent of the lobe 
volume. Radio imaging suggests that the relativistic component is 



clumpy and filamentary, suggesting that mixing has occurred at the 
level of clouds rather than on microscopic scales. 

This evidence of mixing, and the lack of clear cavities raises 
the question of why such a process should be observed in AWM 4 
but not in other systems. Any such discussion is necessarily specu- 
lative. However, one important factor is the age of the AGN out- 
burst. We estimate the synchrotron timescale of the radio lobes 
to be ^--^170 Myr, compared to t he few tens of Myr considered 
typical for FR-I radio galaxies jMack et"al] IT998.) . Dynamical 
age estimates for sources with detec ted ICM cavities are similar 
jPunn et al.ll2005lJ Birzan et al.ll2008h . These estimates may favour 
younger, more powerful radio galaxies whose lobes radiate more 
power at high frequencies and which can excavate larger cavities in 
denser environments. However, even our shortest synchrotron age 
estimate, 66 Myr, is long in comparison, and an age of 170 Myr 
would be fairly unusual. It may be that over such timescales, en- 
ergy losses from the relativistic plasma are sufficient for the con- 
finement of the plasma to weaken, allowing mixing to begin. This 
seems most likely to occur if entrainment slows the jet velocity by 
a large factor. 

The lack of clear cavities is only apparent in AWM 4 be- 
cause of the availability of both a deep X-ray observation and 
low-frequency radio data. Most studies of cavities have used ra- 
dio observations at frequencies > I GHz, and have identified many 
instances of ghost cavities. X-ray strcutures with no radio coun- 
terpart. Using lower radio frequencies extends the timescale over 
which the radio component remains visible, and this may be al- 
lowing us to observe systems which have aged to the point where 
the cavities break up. However, deep X-ray observations are also 
needed, so that even at large radii, where the expected surface 
brightness contrast is small, we can be certain that the cavity is 
weaker than expected. 

There is some evidence of similar "missing" cavities in other 
systems. One example is HCG 62, in which low-frequency ra- 
dio observations rev eal emission exte nding beyond the cavities de- 
tected in the X-ray jGitti et alj|2O10l) . This suggests the presence 
of a set of old, outer radio lobes for which no cavities are detected. 
There is also the possible 'ancient' bubble in the Perseus cluster 
jPunn et al. I l200a> which is observed only as a temperature struc- 
ture, with no known radio or X-ray surface brightness counterpart. 
The age of this bubble is estimated as 100 Myr, comparable to that 
of the AWM 4 radio lobes. However, both these cases differ from 
AWM 4 in that they appear to be related to old AGN outbursts, 
rather than to ongoing activity. Further investigation is needed to re- 
solve these issues, particularly since it seems likely that an increas- 
ing number of similar cases will be observed as low-frequency ra- 
dio observations become more common. 

The question of whether mixing occurs only in the lobes or be- 
gins through entrainment of unheated material in the jets remains 
unresolved. 4C-I-24.36 has a classic FR-I morphology, with nar- 
row (unresolved) jets in the central few kpc, broadening rapidly 
to a width of several kpc over most of their length, with relatively 
diffuse lobes and no hot spots. Such jet broadening is thought to 
occur because of entrainment o f gas within t he jets, leading to a 
rapid decline of jet velocity (e.g. lWorralll2009l) . Slowing of the jets 
may also lead to instabilities which could explain the "wiggles" in 
the jets, as discussed by GVM08. Our estimated values of k could 
be explained if external gas has been entrained and heated by the 
jets. The large uncertainties mean that we cannot be certain that 
entrainment is required, and entrained heated gas is likely not de- 
tectable without very deep hard X-ray observations. Abundance 
mapping does provide evidence of gas motions along the jet axis 
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( IO'Sullivanetai]|2010t) , but it is unknown whetlier tiie enriched 
material is inside the jets or drawn out alongside them by buoy- 
antly rising lobes iChurazov et all 20 00*). 

The jets and lobes of cluster central radio galaxies are com- 
monl y found to be out of pressure equilibr i um with the ICM 
(e.g., bunn & Fabiaij|2004l ; lOunn et aljr2005l : ISirzan et alj|2008l : 
ICroston et al.ll2008h . The range of pressure differences is extreme, 
with pressure ratios of one to several thousand, and the cause of 
these differences is currently the subject of debate. This is gener- 
ally cast in terms of the k parameter, as higher apparent pressure 
imbalances imply a need for a greater ratio of total energy to mea- 
sured energy in the electron population. 

Examination of cavities and radio bubbles in the Perseus and 
Centaurus cluster shows that their k values increase wit h the radius 
of the cavity from the cluster core, and therefore its age jPunn et al] 
l2005h . and this co rrelation is see n in a small sample of radio sources 
in galaxy groups jCroston et a l. 2008). A correlation is also found 
between k and the synchrotron age of cluster radio ga laxies, deter- 
mined from the break frequency of their radio spectra jSirzan et al.l 
I2OO8I) . These findings can be interpreted as supporting spectral ag- 
ing as the cause of the imbalance, with a larger fraction of the par- 
ticle population of older sources falling to lower energies and radi- 
ating at lower frequencies. 

Alternatively, entrainment of additional non-radia ting articles 
within the jets may provide the additional pressure. iDunn et aU 
ilOOf} ) argue M87 and NGC 1275 have electron-positron jets on 
small scales near the AGN, but lobes with high k values which re- 
quire the presence of protons (or more exotic jet models), which 
must therefore be acquired through entrainment. ICroston et al.l 
( I2OO8 ) show that the apparent pressure imbalance is linked to mor- 
phology, with "plumed" FR-I sources (in which much of the lobe 
is beyond the end of the collimated jets) have greater imbalances. 
As "plumed" sources have a greater surface area in contact with 
the ICM, they interpret the difference in pressure as arising from 
a greater degree of entrai nme nt. This result is strengthened by the 
choice of ICroston et alj ^20081) to assume a minimum particle en- 
ergy of 7min=10, thereby including an estimate of the additional 
energy available from low-energy particles and hopefully reducing 
the influence of spectral aging. On the other hand, the character- 
isation of structures in the X-ray observations of the systems is 
hampered by the low X-ray surface brightness of the systems, and 
estimation of the filling factor from the X-ray data is not possible. 

Our observations of AWM 4 are well suited to a study of the 
pressure imbalance and particle content of the radio lobes. We have 
high-quality radio data extending to low frequencies which allow 
us to place relatively strong constraints on the synchrotron pressure 
and age, and make measurements in multiple regions. We are able 
to place some constraints on the filling factors of the lobes, for a 
source where 7^ 1. Our source is also old and therefore near- 
certain to be in equilibrium with its environment. 

It is clear that the choice of energy range over which we as- 
sume the particle population to follow a power law distribution has 
a dramatic effect on our estimates of pressure and k. Taking the 
standard range of 10 MHz- 100 GHz, we find low synchrotron pres- 
sures, high values of k, and therefore a requirement for entrained 
material. The associated radiative ages of the lobes are approxi- 
mately twice the longest dynamical timescales. It therefore seems 
likely that the energy of the electron population is underestimated, 
and the minimum energy cutoff at 10 MHz is too high. 

Using the revised equipartition conditions and 7„i„=100 re- 
duces the pressure imbalance and the required value of k consider- 
ably, but does not reach equilibrium. Our best estimate of the syn- 



chrotron timescale under these conditions (~170 Myr, see TableO 
is rather longer than the buoyant and refill timescales. However, 
the low filling factor of the lobes may bring these timescales into 
agreement. If we lower the minimum energy to ^^^^=10, pressure 
equilibrium in the lobes can be achieved, with k values consistent 
with minimal (or no) entrainment of external material. However, 
this may imply a radiative age for the lobes short enough that purely 
buoyant expansion of the source is ruled out and supersonic expan- 
sion over some part of its lifetime is required. This is not physically 
unreasonable, and may be desirable for energetic reasons. 

We conclude from these results that a large fraction of the ap- 
parent pressure imbalance is due to spectral aging and that inclu- 
sion of low-energy particles by extension to lower values of ^min 
produces more realistic results. We are then left with the ques- 
tion of whether we consider the shorter timescales associated with 
7mm=10 to be likely. 

The relatively flat pressure profile of the jets and lobes (cal- 
culated for 7mi„ = 100) suggests that entrainment from the ICM is 
not effective over most of the length of the jet. In this case, if en- 
trained material does contribute to k it must enter the jet close to its 
source, on scales too small for us to examine. This could suggest 
that entrainment and heating of external material is only effective 
while the jet is collimated and has a high velocity, in its first few 
kpc in the core of NGC 605 1 . This makes supersonic expansion of 
the jets on scales of tens of kiloparsecs less likely, and argues for 
outburst timescales comparable to the buoyant rise time and higher 
values of k. This is certainly the more conservative solution. 

We can also consider the thermal state of the ICM. Our esti- 
mate of the total enthalpy of the radio lobes suggests that it is in- 
sufficient to reheat a large cool core comparable to that in MKW 4. 
However, a smaller core would require less energy, and additional 
energy may be available from shocks or other forms of heating. Us- 
ing our best estimate of the outburst timescale, ~170 Myr, the me- 
chanical energy of the jets is comparable to or slightly less than the 
rate of energy loss through radiative cooling. This could indicate 
that heating was greater at earlier times, or that a shorter outburst 
timescale is required. For the ~66 Myr timescale associated with 
'y^i„=10, the observed filling factors would be more than sufficient 
to balance radiative losses from the ICM, and the supersonically 
expanding jets would certainly drive shocks. 

As discussed at the beginning of this Section, under either sce- 
nario the AGN outburst has a relatively long lifespan. The pres- 
ence of a corona in NGC 6051 may provide some explanation. This 
small volume of cool gas appears capable of fuelling the AGN out- 
burst, replenished by stellar mass losses. As the conduction of heat 
from the surrounding ICM is suppressed, the corona is likely mag- 
netically separate from its environment. This implies a breakdown 
of the ICM-AGN feedback loop, since gas cooling from the ICM 
cannot reach the central engine, and the AGN jets do not heat the 
corona significantly. If this is the case, the AGN outburst is not 
self-limiting and could continue indefinitely. 

The origin of the corona is unclear. ISunl ( l2009l) suggests that 
powerful radio AGN may destroy large cool cores in galaxy group 
scale systems, leaving only the corona in the central galaxy. The 
size of the corona may then be determined by the radius at which 
the radio jets broaden and begin to interact with the ICM, or might 
itself be the cause of this broadening, with the jets losing colli- 
mation as they cross the strong density gradient (or possibly mag- 
netic field) at its boundary. Associations between jet broadening 
and changes in density have been found in other corona-type sys- 
tems (e.g.. Sun et al. 2005a;b). 

In the first case the origin of the magnetic separation be- 
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tween corona and ICM is unclear unless conduction was suppressed 
throughout the pre-existing cool core. Strong suppression of con- 
duction is expected in c ool cores owing to the heat-flux-driven 
buoyancy instability (HBl JParrish & Ouataert 2008). This suggests 
that where magnetic fields are aligned radially along the tempera- 
ture gradient of a cool core, convection will lead to gas motions 
which cause the magnetic field to be realigned tangentially, reduc- 
ing conduction by large factors. It is unclear what the timescale 
for the saturation of the instability would be in a relatively poor, 
low-temperature system such as AW M 4 (the timescale fo r ~5 keV 
clusters is estimated to be ~120 Mvr. lParrish et al.ll2009l) . but HBI 
provides a potentially feasible mechanism. 

The corona could also represent the original galaxy halo of 
NGC 6051, compressed and/or stripped by the surrounding cluster 
gas, with the galactic magnetic field maintaining separation from 
the ICM even before the AGN outburst. The temperature of the 
material outside the corona prior to the outburst cannot be known in 
this scenario. In both cases, the separation between corona and ICM 
suggests that metals released from stars within the corona cannot 
easily diffuse outwards, and that any enrichment must be driven by 
losses from the stellar population outside the corona. 

There is also the possibility that the corona has been stable 
over a timescale considerably longer than the current AGN out- 
burst, and that AWM 4 has been relatively isothermal for a long 
period. If the corona is a stable feature, it is difficult to predict 
the AGN duty cycle, since ICM cooling would not be the driv- 
ing force in triggering activity. The merger of another galaxy with 
NGC 6051 might disrupt the corona, but as NGC 6051 is much 
larger than any other cluster member galaxy and appears undis- 
turbed, it seems unlikely that a merger has occurred for at least 1 
Gyr. For radiative cooling in the ICM to be balanced by energy in- 
jection from the AGN without a feedback mechanism linking the 
two would be highly coincidental. As the roughly isothermal tem- 
peratures demonstrate that the AGN is providing at least enough 
energy to prevent cooling, we would expect that its energy output 
over the long-term would exceed radiative losses, and that no cool 
core has been able to develop. The relatively steep rise in gas frac- 
tion at ~75 kpc radius reported by OS05 could be an indication that 
strong heating has moved ICM gas out of the cluster core. 

A final consideration is the motion of NGC 605 1 within the 
cluster. Both the X-ray morphology and the bending of the radio 
jets suggest that the galaxy is moving south relative to the ICM, 
and the difference in size between the east and west jets suggests 
that there may be an eastward component to the motion. The galaxy 

S iopulation of the cluster is primarily aligned on a north-south axis 
Koranvi & Gelled l2002h as is the ICM, so motion in this direc- 
tion is perhaps unsurprising. Movement of NGC 605 1 relative to 
the ICM may act to spread the energy released by the AGN more 
evenly through the cluster core, and aid in mixing enriched gas out- 
ward from the galaxy. 



8 CONCLUSIONS 

We have used a deep, ~75 ks Chandra observation of the poor 
cluster AWM 4 to examine its structure and properties, and the re- 
lationship between the central radio galaxy and the ICM. Previ- 
ous studies of AWM 4 found the cluster to have a number of un- 
usual and conflicting features. GMRT observations showed that its 
dominant galaxy hosts an old, active FR-I radio galaxy, but XMM- 
Newton found no evidence of cooling in the cluster core to fuel 
this AGN. Conversely, heating a cool core to produce the approx- 



imately isothermal ~2.6 keV ICM observed required more energy 
than was estimated to be available from the radio source. Our analy- 
sis provides solutions to some of these problems, as well as insights 
which may be applicable to other clusters and cluster central radio 
sources. Our results can be summarised as follows: 

(i) The Chandra observation reveals a small cool core located at 
the centre of NGC 605 1 and coincident with t he core radio source. 
This meets the criteria for a galactic corona ( ISunetal.|[2007l : ls^ 
I2OO9I) . It is compact (radius ~ 1-2 kpc), significantly cooler than the 
surrounding cluster halo (A;_Br=1.0lo '^g keV compared to ~2.6 
keV for the ICM), and has a short cooling time (ISllg^** Myr) and 
moderate luminosity (Lx, 0.5-2 = 1-76 x 10*" erg s^^). Heat con- 
duction at the Spitzer rate would be sufficient to heat the core to the 
temperature of the surrounding ICM in 10-20 Myr. This suggests 
that conduction is strongly suppressed. Similarly, a few percent of 
the mechanical energy of the radio jets would be sufficient to have 
heated the corona over the lifetime of the AGN outburst, and we 
conclude that any interaction between the jets and corona must be 
extremely inefficient. VLA 4.9 GHz radio maps do not resolve the 
jets inside the corona, suggesting that they are collimated and nar- 
row, only broadening at its outer edge. 

(ii) We estimate the mass deposition rate through radiative cool- 
ing of the corona gas to be Afcoo;=0.067 A/eyr"^. This would be 
sufficient to power the AGN, requiring an efficiency in convert- 
ing the cool gas to energy of only 0.1 per cent. Mass loss from 
stars within the corona appears sufficient to approximately balance 
cooling losses from the corona. Direct accretion from the 1 keV 
gas at the Bondi rate could fuel the AGN, though the accretion 
rate is rather uncertain owing to the large extrapolation in radius 
required. Magnetic separation of the corona from the ICM would 
prevent gas cooling from the ICM reaching the central engine, and 
the AGN jets do not significantly heat the corona. These factors ap- 
pear to preclude a feedback relationship between AGN and ICM. 
However, the corona itself is capable of fuelling the AGN for long 
periods, and this may explain the unusually long outburst timescale 
estimated for the radio source. It may also explain the lack of a 
large cool core in the system, since without a feedback relation- 
ship, AGN heating seems likely to have exceeded cooling over the 
recent history of the cluster. 

(iii) Imaging shows the gaseous halo of AWM 4 to be generally 
relaxed, in agreement with previous observations. There are weak 
indications of structures associated with the radio source, includ- 
ing a broad bay-like structure around the west jet and lobe. How- 
ever, the only statistically significant surface brightness feature is 
a decrement near the centre of the east lobe. There is no evidence 
of spectrally hard emission associated with the lobes, and the ex- 
pected level of inverse Compton emission is below our detection 
threshold. If the lobes contained only relativistic plasma, we would 
expect to detect the cavities with high statistical significance. We 
interpret these results as indicating that the lobes are only partially 
filled. This is supported by the clumpy, filamentary appearance of 
the lobes in radio images. Based on the surface brightness, we can 
place limits on the fraction of ICM plasma in the lobes. Assum- 
ing the remaining volume is occupied by radio-emitting relativistic 
plasma, we find filling factors for this component of = 0.24 and 
= 0.21 for the east and west lobes respectively, with 3a upper 
limits of(j}< 0.43 and (t> < 0.76. 

(iv) We measure the pressure profile of the ICM, and compare 
this with minimum energy pressure estimates for the jets and lobes 
of the radio source. Under the most conservative conditions, as- 
suming contributions from particles emitting between 10 MHz and 
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100 GHz, we find a strong pressure imbalance between the lobes 
and their environment, with the lobes apparently underpressured 
by a factor ~160. However, these conditions imply an age for the 
source which is considerably longer than the timescale for the lobes 
to buoyantly rise to their current position. Estimates which include 
less energetic particles greatly reduce the pressure imbalance, to a 
factor ~15 for 7^i„=100, or to a factor ~4 for ^^i„=lO, with the 
latter case consistent with pressure equilibrium within uncertain- 
ties. The radiative age estimated for 7„i„ = 100 is roughly consis- 
tent with the buoyant timescale of the lobes, while the age estimated 
for 7,„i,j=10 is significantly shorter and could imply a supersonic 
expansion phase for the jets. We consider the effects of bending in 
the jets on our pressure estimates, but find that they are unlikely to 
have a significant impact on our results. 

(v) From these measurements and the filling factor estimates de- 
scribed above, we estimate the required ratio of energy in non- 
radiating particles to the observed energy in electrons, k, for each 
lobe. Assuming ymin ~ 100, we estimate fc=37.5 and 24.8 for 
the east and west lobes, with 3a upper limits of fc < 741.6 and 
k < 517.0 respectively. For "fmin = 10, we estimate fc=10.4 and 
A;=6.8 for the east and west lobes, with large uncertainties consis- 
tent with k=l (an electron-proton plasma) or fc=0 (an electron- 
positron plasma), or with the ymin = 100 values. This indicates 
that in principle the apparent pressure imbalance in the lobes can 
be resolved by the inclusion of these lower energy relativistic parti- 
cles. Alternatively, entrainment and heating of thermal plasma (ei- 
ther from stars or the ICM) could provide the necessary additional 
pressure. However, such entrained material must have a low den- 
sity and very high temperature, rendering it undetectable with the 
current data. 

(vi) We estimate the enthalpy of the radio lobes and find that for 
the filling factors estimated above it is ~ 0.3 — 7.0 x 10^** erg. This 
is lower than the estimated energy required to reheat a large cool 
core such as that seen MKW 4, a cluster of similar temperature and 
galaxy population. The mechanical power output of the jet depends 
on the timescale of the outburst; for our best estimate of ~170 Myr 
the jet power is ~ 0.6 — 1.3 x 10*'^ erg s~^. This is lower than or 
comparable to the bolometric X-ray luminosity of the ICM within 
the radius of the lobes, ~ 1.3 x lO'*'^ erg s~^, suggesting that in the 
absence of other forms of heating, the energy available from the 
radio lobes is at best just sufficient to balance cooling in the ICM. 

However, the relative isothermality of the ICM and lack of any 
significant cooling region outside the corona strongly suggests that 
the AGN has provided enough energy to balance or exceed cooling 
losses in the past. This could be achieved through additional heating 
mechanisms (e.g., weak shocks, sound waves, cosmic rays), larger 
filling factors, or a shorter outburst timescale. If we instead assume 
the shorter radiative timescale (~66 Myr) estimated from the spec- 
tra of the lobes and assuming 7mi„=10, the jet mechanical power 
is 2.6x10*"^ erg s~^, in excess of the cooling rate. This timescale 
would also require supersonic expansion of the jets, providing ad- 
ditional heating through weak shocks. 
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